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Abstract 

For the quintom models with arbitrary potential V = V(<p,cr), the 
asymptotic value of equation of state parameter u is obtained by a new 
method. In this method, ui of stable attractors are calculated by using the 
ratio dlnU/dlno in asymptotic region. All the known results, have been 
obtained by other methods, are reproduced by this method as specific 
examples. 



1 Introduction 

The recent observations on type la supernova Tj, large scale structure [2] and 
cosmic microwave background radiation [3] , support the claim that our universe 
has accelerating expansion and it is composed of dark energy (73%), dark matter 
(23%) and baryonic matter (4%). 

The simplest candidate for dark energy is a cosmological constant A of order 
(10 _3 ev) 4 , but it suffers from conceptual problems such as fine-tuning and coin- 
cidence problems [I]. As the alternative to cosmological constant, the dynamical 
models have been introduced. 

One of the important parameters of dark energy is the equation of state pa- 
rameter oj = p/p, which must satisfy lu < —1/3 if one interested in accelerating 
universe. The time dependence of w, i.e. cj(z), is usually used to explain the dy- 
namical nature of our universe. Some astrophysical data seems to slightly favor 
an evolving dark energy and shows a recent ui = — 1, the so-called phantom- 
divide- line, crossing [5]. 

Many of dynamical dark energy models are based on scalar fields. The 
simplest of these models is the quintessence model which consists of one normal 
scalar field <p [BJ. In all quintessence models, always u> > — 1. The other kind 
of scalar field model is phantom model. In this model, there exists a phantom 
scalar field a which its kinetic term appears with minus sign. In phantom model, 
u) always satisfies uj < — 1 |7J. None of these two models can explain the u> = — 1 
crossing and at least two scalar fields are needed, in models known as hybrid 
models, for occurrence of this transition. One of these models is the quintom 
model which consists of one quintessence and one phantom fields [5]. One can 
show that in all quintom models with slowly-varying potentials, the transition 
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from oj > — 1 to cj < — lis always possible [9 . Related to oj = — 1 crossing, it 
must be added that the occurrence of this transition is also possible in a model 
with only one scalar field, but this field must be coupled to background matter 
field [ID]. 

One of the interesting topics in dark energy problem is the studying of the 
late-time, or asymptotic, behavior of physical quantities. In this connection, the 
investigation of attractors of the various dark energy models is one of the main 
tools [H]. Among the attractors of a dynamical system, the stable attractors 
determine the late-time behavior of various quantities, including the equation of 
state parameter oj. Despite its usefulness, the attractors' studies are restricted 
to specific potentials. This is because the choosing of the set of convenient 
variables, which results in an autonomous set of equations, is a highly nontrivial 
task and depends on the functional form of the considered potential. Therefore 
seeking other methods for studying the late-time behavior of, for instance, oj is 
an important task. 

In ref.[12j. the asymptotic behavior of oj a of phantom models has been in- 
vestigated for potentials which V(a) — > oo asymptotically. There, it has been 
shown that all such models are divided into three classes, characterized by 
oj a — > — 1, oj a — > ojq < — 1 and oj g —* — oo. In ref.[13j. the evolution of a special 
quintom model, the so-called hessence model, has been discussed in oj~ oj' plane. 
The hessence model is a quintom model with potential V((p, u) — V((p 2 — a 2 ), (p 
and a denote the quintessence and phantom fields, respectively. In this model, 
there is, effectively, only one dynamical field <p = \/ tp 2 — a 2 ■ The late time 
behavior of this model has been classified in oj — oj' plane in 13J and it has been 
shown that there exist four distinct regions in this plane. Both of the above 
mentioned studies are restricted to cases with only one dynamical dark energy 
field. 

In this paper we want to study the asymptotic behavior of equation of state 
parameter oj for quintom models with arbitrary potential. Since it is a two- 
field-component model, there exists new features which do not appear in single- 
component models of [12] and [E]. It is must be added that the attractors of 
quintom models with two classes of potentials has been recently studied in [14j . 
The potentials are restricted to cases in which the scalar fields do not diverge, 
and the potentials which behave asymptotically as an exponential potential. 
As we will see, there are too many examples that do not belong to the above 
mentioned cases, (to be specific, the examples B, C, D, F and G of section 2), 
and one must seek another method, as we do, to obtain the asymptotic behavior 
of the system. 

The paper is organized as follows. In section two, we briefly review the 
quintom model and derive a relation which describes the dynamical evolution of 
the quintom potential. This equation is one that we use to obtain the asymptotic 
behavior of equation of state parameter oj of an arbitrary quintom model. It is 
the same relation which exists in hessence model |13j , and reduces to the known 
relations in quintessence [15j and phantom [16] models in appropriate limits. 
In section three, we use this relation to obtain the asymptotic value of oj for 
several examples. The examples are divided to two classes: the potentials which 
have no quintessence-phantom interaction, i.e. V(ip, a) = V±(ip) + V2(cr), and 
the potentials containing the interaction term V[ n t(y>, a). In all examples, it is 
seen that our results coincide with those obtained by other methods, e.g. the 
attractors' study or the numerical study of fields' equations. 
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We use the units c = h = 8itG = 1 throughout the paper. 



2 The evolution equation of quintom potential 

we consider a spatially flat Friedman- Robertson- Walker space-time in comoving 
coordinates (t, x, y, z) 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (1) 

where a(t) is the scale factor. The quintom dark energy consists of the quintessence 
field (p and the phantom field a, with Lagrangian 



Ac = \d^d^ - \d^a - V(<p, a). (2) 



The energy density pde and pressure pdc of the homogenous quintom dark energy 
are 

Pdc =~<P 2 - ^o- 2 + V(<p,a), 

Pdc =~¥> 2 - ^o- 2 - V(ip, a), (3) 
and the evolution equations of the fields are 

+ 3*0 + ^2=0, (4) 

and 

t + 3H*-^p^=0. (5) 
do- 
Here H(t) = a(t)/a(t) is the Hubble parameter and 'dot' denotes the time 
derivative. The Friedman equations are 

3H 2 = Pdc , (6) 

and 

- 2H = pdc +Pdc- (7) 

Note that the equations (J4j)- (JTJ) are not independent and eq.Q is obtained from 
eqs.Q-©. The equation of state parameter u> = Pde/ Pde is therefore 

ip*-& 2 -2V(<p,a) 

Because of the field equations Q and ((SJ), the quintessence field ip rolls down 
the potential to reach a minimum of V, while the phantom field a falls up the 
potential V and settles in a maximum of V. So the late-time values of (ip, a) 
are the saddle point of potential V, we denote it by (<£>*, cr*). We will use this 
fact in our analysis. 

To obtain the evolution equation of potential V, we first define the function 

1 — U) 2V 
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if cj 7^ — 1. Differentiating In a; with respect to In a results in 




(10) 



Using eqs.(g]) and © and noting that V = (dV/d(p)(p + (dV/da)&, eq.([TU| can 
be written as 

dlnx 1 f ... 2 V 



dlna H I 1+wU 

or 



(11) 



-7-w+*H*£)- (12) 

This is a general equation which can be used for ip = case (which the quin- 
tom model becomes phantom model), a — case (where the model is the 
quintessence model), and when V(ip,a) = V{4> = \/ f 2 — & 2 ), which the quin- 
tom model becomes hessence model. In all above mentioned cases, the potential 
is a function of only one field variable, e.g. the field 9, and therefore eq. (fT2|) can 
be reduced to an equation for the field-derivative of potential. This is because 
V — (dV/d9)6. But in the quintom model, it is not the case. 

It is more suitable to write the left-hand-side of eq. (fT2^) in terms of derivative 
with respect to In a. The result is 

r 3(1+0,) (l+l-), (13) 



V " ' \ 6 
in which 

f - = U d4) 

dm a H 

We are going to use eq. ()13|) to determine the late-time behavior of to for quintom 
models, i.e. to obtain the value uiq for stable attractors of quintom models. To 
do so, we must first evaluate the ratio x' fx for attractor solutions. The first 
step in studying the attractors is to introduce a set of convenient dimcnsionlcss 
variables, which in study of quintom models, three of these variables are 



y=-nTT T - ( 15 ) 



The other variables depend on the functional form of the potential, but the 
above three variables always exist. In terms of these variables, the ratio x'/x 
becomes 

x[_ _ 2iu' _ 2 [x v y 2 x' v - x a y 2 x' a + (x\ - x%)yy'] 
x 1 — uj 2 y 2 {x 2 — x 2 T ) 

For attractor solutions, one has x' — x' a = y' = 0, so x'/x — > 0. Of course 
one must be more careful for situation in which y 2 (x 2 P — x 2 ) = 0. In quintom 
models, the fields tp and a satisfy different equations of motion, i.e. ip moves 
toward the minimum of potential while a goes to the maximum of V. Therefore 
we expect that in stable attractor solutions, y 2 ~ V(ip*,<r*) ^ and usually 
xz, 7^ x\. In next section, we see that these conditions are fulfilled in all known 
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examples, except the example F which we discuss it separately. So for stable 
attractors, the equation of state parameter oj tends to ljq, obtained from 



1 V 

-o = -l-3 Y - (17) 



3 Examples 



In this section we apply the main result (I17p to situations which have been 
studied by other methods. These examples are divided to two classes: The 
quintom models without quintessence-phantom interaction and the potentials 
contain this interaction term. We discuss these two kinds of potential separately. 

3.1 The potentials V((p, a) = Vi(<p) + V 2 (a) 

For potentials of the form V(<p,a) — Vx{<p) + V 2 (a), the asymptotic configura- 
tions of the fields p and a are (<p*, <x*), where Vi(ip) becomes minimum at ip* 
and V%(a) attains its maximum at a*. In following we consider four examples. 

A. V=V (po e- x ^ + V ao e- x 2° 

The attractors of this potential have been studied in [T7j. As the exponen- 
tial function has no extremum points, so we can only use the fact that ip always 
rolls down the potential and a falls up it. So <p* — > oo and a* — > — oo, from 
which V(p, a) — > V 2 {cr). The ratio V'/V then becomes 

v ^vi {dv 2 /day , 

V^V 2 = — v 2 = - X2<T - (18) 

Eq. (fT7|) then results in 

u luq = -1 + ^A 2 cr'. (19) 

This is exactly the result obtained in [T7]. The only stable attractor of this 
potential is characterized by x v = 0, x a — ~\ 2 /V6 and y 2 = 1 + A|/6. So 



V 2 _ _ A 2 
xl - x 2 a + y 2 3 



1 - -7T, ( 20 ) 



which is equal to p^|) ( note that x a — a'/y/6 = —\ 2 /V6, results in a' = —\ 2} 
sow = -l- Xl/3 = -1 + (l/3)A 2 cr')- 

B. V = V Vo p a + V ao a a (a>0) 

The time- variation of u> of this potential has been studied for two special cases 
a = 2 and a = 1.8 in 18J. V\{p) is minimum at <p* = 0, but V 2 {a) has no 
maximum, although it is clear that |er*| — > oo. So V(ip,a) — > V 2 {a) and 



V Vi _ 2a' 



— > 0, (21) 



V V 2 a 
from which eq. ()17j) results in 

LU^LU = -l. (22) 
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This is the same result obtained for special cases in [TH]. It can be shown that 
eg. (|2"2"|) is also valid for odd a's where ip* — * —oo. 



C. V = V Vo e- x ^ +V ao e- x ^ 

Vi((p) has no minimum, but it is obvious that \(p*\ — > oo. V2(cr) is maximum at 
a* = 0. So again V((p, a) — > V~2{<j), from which 

y - y = -2A 2 aa'| CT=0 = 0. (23) 
This results in luq = — 1 which is consistent with numerical calculation of [17] . 

D. V = V Q e Xu " (a > 0) 

Now we consider a phantom model with the above potential. In this case 
1(7*1 — > oo. Therefore 



(24) 

I oo, a > 1 

from which 

f-1, a<l 

(25) 

In above we use the fact that in a=even cases, we have a* — > ±oo. For a* — > 
+oo, since a' > 0, one has <j a ~ 1 a' — * +oo. For a* — * — oo, wc know that 
a' < and therefore again cr^V — > +oo. This general result agrees with the 
numerical calculations of [12] which have been done for special cases a = 0.5, 1, 
and 2 with A = 1. It is clear that eq. (f2"5j) can be also used as the late-time 
behavior of lo of a quintom model with the potential V — V Vo e Xlip + V CTo e A2Cr . 

3.2 The interacting potentials 

We now consider the potentials of the form V(ip, a) — Vi(( / 9) + V r 2(cr) + Vi n t.( ( ^, c), 
through three following examples. 

E. V = Voe~ v/ ^^ mip+ "' T ) 

The attractors of this quintom potential have been discussed in [19] . This 
potential has no extremum point, but it is clear that ip* — > oo and a* — > — oo. 
The ratio V'/V is 

V 1 fdV . dV A r-, . , 
y=y{^'+^')=-^(W + na'), (26) 




so 



Wo = — 1 + — [pup + na ) 
o 

= — 1 + 2(mx v + nx a ). (27) 
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In [13] , it has been found two stable attractors, denoted there by P and T solu- 



tions. In P-case, the coordinates are : 



m, x 



-n and y = 1- 



^0, 



which result in wp = — 1 + 2 (to 2 — n 2 ). This is the same as our result (j2"T|) . In 
T-case, it has been obtained x v — to/(2(to 2 — n 2 )), x a ~ —n/(2(m 2 — n 2 )) and 



y 2 = 1/(4(to 2 — n 2 )), with wt = 0. This solution is also obtained from (|27|) . It 
is interesting that the oj of other solutions of [19] are not reproduced by using 
(f2"T)) . It is because these solutions are unstable attractors and therefore do not 
describe the late-time behavior of the system. 



F. V 



±m 2 <p 2 + \g 2 ^ 2 o 2 + (M 2 - Aa 2 ) 2 /(4A) 



This potential has been investigated in [2]. In this case, it can be easily seen 
that the potential has a unique saddle point at (</?*, <r*) = (0,0). Therefore 



V 



1 fdV 



V \d<p 



dV , 
da 



0. 



(28) 



(v=,<t) = (0,0) 



which results in ujq = — 1 . This is the same value obtained in [14) , denoted there 
by de Sitter solution. Note that here x v — x a — and y = 1. So eq. lflUI leads 
to 

— , (29) 
x y 

which again becomes zero at y' = and y = 1. 

G. V = Vi + V 2 + \VV1V2, where V x = V Vo e- x ^ and V 2 - V^- e ^ 2rT 



This is the potential whose attractors have been obtained in 20J . This po- 
tential has no cxtremum point, but it is clear that ip* — > 00 and a* — > — 00. In 
this case 



V 



1 

V 



X 



X 



(V r i,V 2 )--(0,oQ) 



A 2 a'. 



So 



w 



-1 



1 



A 2 a'. 



(30) 



(31) 



The only stable attractor of this potential is the second solution of Table 1 of 
[2"0] with coordinates x v — 0, x a = — A2/V6 and y 2 = 1 + A 2 /6. These values 
result in w = — 1 — A^/3 = — 1 + A2<t'/3, which is same as eq. (|3~Tj) . 
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